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We study the elasticity of cross-linked networks of thermally fluctuating stiff polymers. As com- 
pared to their purely mechanical counterparts, it is shown that these thermal networks have a 
qualitatively different elastic response. By accounting for the entropic origin of the single-polymer 
elasticity, the networks acquire a strong susceptibility to polydispersity and structural randomness 
that is completely absent in athermal models. In extensive numerical studies we systematically vary 
the architecture of the networks and identify a wealth of phenomena that clearly show the strong 
dependence of the emergent macroscopic moduli on the underlying mesoscopic network structure. 
In particular, we highlight the importance of the polymer length that to a large extent controls the 
elastic response of the network, surprisingly, even in parameter regions where it does not enter the 
macroscopic moduli explicitly. Understanding these subtle effects is only possible by going beyond 
the conventional approach that considers the response of typical polymer segments only. Instead, we 
propose to describe the elasticity in terms of a typical polymer filament and the spatial distribution 
of cross-links along its backbone. We provide theoretical scaling arguments to relate the observed 
macroscopic elasticity to the physical mechanisms on the microscopic and the mesoscopic scale. 

PACS numbers: 87.16.Ka, 62.20.Dc, 82.35.Pq 



I. INTRODUCTION 

Classical elasticity is a continuum theory that deals 
with the large scale deformation properties of solid sys- 
tems. It relates stresses and strains by introducing a 
host of phenomenological parameters, e.g. shear and bulk 
modulus for isotropic media, that characterize the elas- 
tic properties on wave lengths large compared with any 
other material length scale Biological systems like 
the cell or sub-cellular organelles are often characterized 
by a highly heterogeneous structure with a multitude of 
hierarchical levels of organization ■ Due to these large 
scale inhomogeneities that may extend up to the scale 
of the system size, the applicability of elasticity theory 
on smaller length scales has to be critically examined. 
In particular, the actual deformations in the system are 
expected to relate to the externally applied stresses in 
a non-trivial way that crucially depends on the specific 
structural details. 

To shed some light on the relevance of structure to 
the effective elasticity this article deals with the calcu- 
lation of elastic constants in networks of semi-flexible 
polymers. In eukaryotic cells these networks assemble to 
form the cytoskeleton that plays a central role in many 
cellular functions such as locomotion, adhesion or cell di- 
vision. From the point of view of structure already a 
one-component isotropic solution of semi-flexible poly- 
mers represents an interesting model-system being stud- 
ied for many years [H, 0, Q . One of the main quantities of 
interest is the plateau value of the shear modulus found 
at intermediate timescales where single polymer bending 
fluctuations are equilibrated, yet center of mass motion 



is negligible. The generally accepted theory for the con- 
centration dependence of the plateau modulus is based 
on the free energy change of confining a polymer to a 
tube [1, H, 0, @| , the diameter of which is a consequence 
of the structural organization of the tubes in the form 
of a random assembly of cylinders [9( . Even though this 
is well known for more than a decade, computer simula- 
tions to study the geometrical as well as elastic properties 
in this fibrous architecture have only recently been real- 
ized [H El. 

Upon the addition of cross- linking agents or other regu- 
lating proteins one can induce structural changes to mod- 
ify the network architecture in many ways [l3l [Til . 
ITa . HH . There have been attempts to describe the phase- 
diagram of these systems [13, EH , the detailed mecha- 
nisms that lead to a particular structure, however, are 
far from being understood. In general, there will be a 
complex interplay of polymer kinetics, thermal fluctua- 
tions and chemical as well as mechanical properties of the 
polymers and the cross-linking agents yielding a partic- 
ular architecture relevant for a given physical situation. 

A complementary approach to describe cross-linked 
networks is to neglect these intricate "dynamic" as- 
pects of the network, and to concentrate on a "static" 
architecture and its effect on the macroscopic elastic- 
ity [H S3, HI, H, HI Hi- In the structural engineering 
community, for example, it is of tantamount importance 
to analyze the architecture of structures made of beams 
or trusses. A common way to take advantage of the 
reduced weight compared to the bulk material without 
suffering from a loss of stiffness is a triangulation of the 
basic cells. This eliminates the soft bending modes of the 
beams and makes it possible to construct huge cantilever 
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bridges like that over the Firth of Forth in Scotland or 
towers like Eiffel's tower in Paris. Since the rigidity of 
these structures is not due to the individual beam but 
to a non-local back-coupling effect induced by the archi- 
tecture of the network, the triangulation is therefore one 
example on how cooperativity among the building blocks 
may be possible. 

To address this question of cooperativity in the context 
of the elasticity of cross-linked stiff polymer networks we 
will concentrate in the following on two generic struc- 
tures, cellular and fibrous networks, that may serve as 
reference systems for the classification of real polymer 
networks. While cellular structures may be characterized 
by the amount of randomness in size and type of their 
unit cells (see Fig|T}>d), fibrous networks have a hierar- 
chical structure, where smaller cells are generated within 
lager cells within even larger cells (FigQJi). This is a 
consequence of the presence of the additional mesoscopic 
scale of the fiber length. As we will see, this length-scale 
is ultimately responsible for the intricate scaling proper- 
ties of the elasticity of fibrous polymer networks. The 
goal of this article is to identify these mechanisms, that 
couple the particular network structure to the properties 
of the individual polymers and effectuate the macroscopic 
elasticity of the system. 

In contrast to the purely mechanical systems relevant 
for engineering applications [H, HH, EH , the systems 
we would like to study are immersed in a thermal envi- 
ronment. This implies that in addition to the usual en- 
thalpic polymer elasticity also entropic effects have to be 
accounted for. We have published a brief account of this 
study recently [23]. It will turn out that by accounting 
for the entropic origin of the single-polymer elasticity, the 
networks acquire a strong susceptibility to polydispersity 
and structural randomness that is completely absent in 
athermal models. 

The article is organized as follows. In Sect HI] we moti- 
vate our modeling approach of thermally fluctuating net- 
works of stiff polymers. This will lead us to the definition 
of effective elastic properties of the "polymer segments" 
that constitute the elementary building blocks of the net- 
work. In Sect IIIII and IIVI these polymer segments are as- 
sembled into cellular and fibrous networks, respectively. 
The macroscopic elastic constants of these structures are 
calculated and related to the particular architectural fea- 
tures. Finally, in Sect [V] we present our main conclusions 
and hint at implications for experiments. 



II. MODEL DEFINITION 

To study the elastic properties of thermally fluctuating 
cross-linked stiff polymer networks we calculate numeri- 
cally the low frequency shear modulus. Assuming a time- 
scale separation between the fast bending fluctuations of 
the single polymer and their very slow center of mass mo- 
tion, we adopt a description of the system in the spirit of 
a Born-Oppenheimer approximation. This neglects en- 
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FIG. 1: Illustration of the different architectures of (a) fibrous 
and (b-d) cellular materials in two dimensions. While (a) 
and (b) are random structures generated by Poisson point 
processes, (c) and (d) are quite regular networks based on 
honeycomb and square lattices, respectively. 



tropic contributions from the "slow variables" , the cross- 
link positions, while assuming the "fast" polymer degrees 
of freedom to be equilibrated at all times. Macroscopic 
quantities will then depend parametrically on the set of 
cross-link variables. A macroscopic shear strain 7 con- 
strains the cross-links at the boundaries, while those in 
the bulk are moving freely to minimize the elastic energy 
E. The shear modulus is defined as its second derivative 
with respect to the shear strain, G = V~ 1 d 2 E m i n /dj 2 , 
where V is the system volume. 

By keeping the positions of the cross-links fixed, the 
energy can be written as a sum 



(1) 



over contributions from individual polymer segments a, 
each of which connects a given pair of cross-links (see 
Figf2|). The single segment energy e depends on the gen- 
eralized "displacement- vector" Sx a , which incorporates 
the degrees of freedom, displacements u and rotations 6, 
of the two cross- links pertaining to the segment. 

In the numerical section we focus on two-dimensional 
systems such that a vector 5x2d — (uq. 9q,ui,9i) has six 
components. Those are in-plane displacements Uo,; and 
z-axis rotations #0,2 > at both ends 0, 1 of the segment with 
length I (see Fig 12]). Note, that the additional variable 
of cross-link rotation sets our system apart from bond- 
bending and related models [26[ where only translational 
degrees of freedom are accounted for. As a consequence 
one also has to account for the presence of torques as the 
conjugate variable to rotations. 



3 



CL 



^ bending i 



segment of length 1 9 I L 



CL . 



mech stretching 



therm stretching 

K 



FIG. 2: Illustration of a polymer segment of length I and its 
connection to the network (dashed lines) at the cross-links 
(CL). The three degrees of freedom at each cross-link are de- 
noted by u and 9, respectively. Identification of the three 
possible modes of deformation and their stiffnesses k±, k 3 
and fen as defined in the text. 



To leading order, in linear elasticity, the single segment 
quantity e is a quadratic function of its coordinates 



e(x) = jx T Kx, 



(2) 



which defines the "stiffness matrix" K (spring constants) 
of the polymer strand. 

In models of classical beams with cross-section radius 
r the matrix elements are well established and relate to 
the two deformation modes of stretching (s) and bending 
(_L), respectively. While the former is characterized by 
the Young's modulus E of the material, the latter de- 
pends on the bending stiffness k = Eirr 4 /4, here taken 
for circular cross-sections. To calculate the bending re- 
sponse, standard Euler-Bernoulli beam theory [l[ is used. 

While we refer to Appendix [Al for a derivation of the 
complete matrix, it turns out that the response of a beam 
of length I is sufficiently characterized already by two 
elements of K, 



k s {l) = An/lr 2 



k±{l) = 3k/1 3 



(3) 



relating to either deformation mode. Due to their small 
aspect ratios r/l <C 1 slender rods are highly anisotropic 
and much softer in bending than in stretching, k±/k s oc 
(r/l) 2 . In this approximation the two deformation modes 
are decoupled such that, for example, pre-stretching does 
not influence the bending response. Therefore, Euler- 
buckling cannot be accounted for. 

Here, we consider thermally fluctuating stiff polymers 
immersed in a heat bath of solvent molecules. In these 
systems, the effects of temperature on the elastic prop- 
erties of the polymer can be quantified by defining the 
persistence length l p as the ratio of bending stiffness to 
thermal energy l p — K/ksT. With this definition we 
have, in addition to the enthalpic stiffness of the classical 
beam, an entropic contribution 



k\\(l) = C« 



lp 



(4) 



to the polymer's stretching compliance that can be cal- 
culated within the wormlike chain model [27], [28| . The 



prefactor ( depends on the specific boundary conditions 
chosen at the ends of the polymer segment. Its value 
can be absorbed in the persistence length, and therefore 
only quantitatively affects the results. To avoid a large 
numerical offset with respect to Eq.([3]), we have chosen 
C — 6, which corresponds to a boundary condition with 
one end clamped [27] . Having two longitudinal deforma- 
tion modes k s and k\\ the effective stretching stiffness is 
equivalent to a serial connection 
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(5) 



Thus, the elastic properties of the polymer segments 
are given by the classical theory of beam bending sup- 
plemented by a generalized stretching stiffness, that also 
includes entropic effects. While the stiffness matrix has 
only been set up for the two-dimensional problem, the 
governing entries in three dimensions will still be the 
same Eqs.© and (gj). 

As one can infer from Eqs.© and at a given tem- 
perature T there are two length scales characterizing the 
material properties of the polymers, the radius r and 
the persistence length l p . Typical biological polymers 
are characterized by a ratio R = l p /r 3> 1. F-actin, 
for example, a key component of the cytoskeleton has 
R = O(10 4 ) (r Ri 5nm , l p w 17/im), while microtubules, 
most important for cell-division and infra-cellular trans- 
port, have an even larger R = O(10 6 ). For specificity, we 
require in the following a constant R — 1.5 • 10 4 , the pre- 
cise value, however, is irrelevant if one is interested only 
in the thermal response where the radius does not enter 
and k s — > oo [45[ . Occasionally, we will perform this limit 
to highlight features that are independent of the mechan- 
ical stretching response. On the other hand, the location 
of the cross-over point, where the mechanical stretching 
becomes relevant, does indeed depend on the choice of 
R. By definition, it determines the relative magnitude of 
the two stretching compliances k s /k\\ ~ R 2 (l/l p ) 3 . 

The dependence of the three force constants k± , k s and 
feii, Eqs.J3]) and (H]), on the ratio of persistence length 
to segment length l p /l is illustrated in FigJ3] One can 
clearly distinguish three regimes, in each of which one of 
the spring constants is by far smaller than the remain- 
ing two. The dashed line corresponds to a hypothetical 
spring where the deformation modes are coupled in series 
k^ 1 = kj 1 + k^r 1 + k^ 1 . If the segment length I was rep- 
resentative for the network under consideration, that is 
the network was characterized by only small polydisper- 
sity, then we would expect the macroscopic modulus to 
be well approximated by the microscopic single segment 
behavior considered here. We will later refer to this be- 
havior as the "affine model" . It will be shown to be valid 
only in regular cellular structures. 

This completes the specification on the microscopic 
scale of the elastic properties of the single polymer seg- 
ments. We now proceed to assemble the segments into 
networks of varying architecture to identify the physical 
principles which determine the elastic response on the 
macroscopic scale. 
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FIG. 3: Dependence of the three spring constants fey , k±_ and 
k s on persistence length l v /l. The dashed line corresponds 
to a hypothetical spring with the three deformation modes 
connected in series. 



To determine the elastic shear modulus, we apply a 
small shear strain of 7 = 0.01 to stay in the regime of 
linear elasticity and use periodic boundary conditions on 
all four sides of the simulation box. The numerical proce- 
dure is performed with the commercially available finite 
element solver MSC.MARC. The results will be comple- 
mented by scaling arguments. 

We will find that in regular cellular architectures, to 
be discussed next, macroscopic elasticity can trivially be 
explained by the microscopic constitutive laws given in 
terms of the stiffness matrix K. In sufficiently random 
cellular systems, however, this picture is changed. The 
macroscopic response takes up nontrivial features that 
cannot be explained by single polymer elasticity. In fi- 
brous architectures, subject of Sect llVl we will find this 
anomalous elasticity again but in more striking form. 



III. CELLULAR ARCHITECTURE 

A cellular structure is most conveniently constructed 
from a Voronoi tessellation of a distribution of points 
which may either be chosen regularly or by some random 
process [29| . With each point we associate a Voronoi cell 
that is defined to enclose that region in space which is 
closer to the given point than to any of its neighbors. 
This procedure is equivalent to the Wigner-Seitz con- 
struction known from solid-state physics. In three dimen- 
sions the elastic elements are defined to be the lines of 
intersection of two neighboring cell walls, while in two di- 
mensions (see FigfTJ they are represented by the cell walls 
themselves. We will call these elastic building blocks of 
the network polymer segments and associate to them the 
material properties, respectively the stiffness matrix K, 
introduced in the preceding section. By its definition, a 
segment spans the distance between two vertices and is 
therefore "end-linked" to the rest of the network. 

Depending on the spatial distribution of Voronoi 



points there will also be a distribution P(l s ) of segment 
lengths l s . Only in regular structures, for example the 
(anisotropic) two-dimensional honeycomb structure, this 
distribution will degenerate into one (or several) delta- 
function peaks. 

The first moment of this distribution, the average seg- 
ment length l s , is naturally the most important quantity 
to describe the geometrical aspects of a cellular struc- 
ture. In d = 2,3 dimensions this "mesh-size" may be 
reparametrized in terms of the density p as 

r.ocp-VW-D, (6) 

where we defined p as the total polymer length per sys- 
tem size. While there are practical reasons to use p as a 
measure for the density in the simulations, in experimen- 
tal work it is sometimes easier to control the monomer 
concentration c. This can be found as rc oc p, where 
the cross-section radius r is assumed proportional to the 
monomer size. 



A. Mechanical Behavior: Beams 

In the engineering literature the cellular structures de- 
fined above are well known as foams and are ubiquitous 
in nature and many areas of technology. Examples range 
from liquid foams and froths well known from drinks 
or household detergents, to plastic and metallic foams 
used for insulation or shock absorption [29|, |3fJ. It is 
well known that naturally occurring foams have to obey 
Plateau's laws to reach an equilibrium state. We do not 
require these laws to hold in the following, since we are 
interested in the dependence of elastic properties on the 
architectural features in general, and not in the specific 
details of the dynamic properties of foams. 

For purely mechanical cellular foams, where thermal 
fluctuations are neglected altogether, the only material 
length scale is the radius r of the cross-section. By iden- 
tifying k/I s as an energy scale, we can use dimensional 
analysis to write the shear modulus G as 

G=J^g(r/T a ), (7) 

where the occurrence of the spatial dimension d high- 
lights the fact that the modulus has units of an energy 
density. In writing this, we have not made explicit the de- 
pendence on the higher moments of the probability distri- 
bution P. As will become clear below, these can be used 
to characterize the randomness of the structure and will 
be considered separately. If one defines force-constants 
at the scale of the average mesh-size 

k ± ~ k/IJ , k s ^ n/l s r 2 , (8) 

the scaling variable can alternatively be written as r/l s ~ 
yk±Jks, and therefore characterizes the relative stiffness 
of the bending to the stretching mode. 
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B. Regular Structures and Affine Models 

Restricting our attention for the moment to regular 
structures, macroscopic elasticity can already be under- 
stood by considering the response of a single cell [3(| 
l3lL I32I ]. In these systems it seems reasonable that lo- 
cal stresses acting on an individual cell are the same as 
those applied on the macroscopic scale. In other words, 
the local deformation 5 of a cell with linear dimension l s 
follows the macroscopic strain 7 in an affine way such 
that it scales as 5 oc jl s . With this assumption the scal- 
ing function can be calculated [3l[ and one generically 
finds for the modulus 




G-^l^iak-'+bk- 1 ). 



(9) 



where the details of the particular structure may enter 
the numbers a and b in an involved way. The important 
conclusion to be drawn is that the deformation modes act 
as if they were springs connected in series. For slender 
beams with r < ! s the bending mode is softer than the 
stretching mode and therefore dominates the modulus - 
mechanical foams are bending dominated. 

While we argue here that the modulus in Eq.© rep- 
resents the generic case, there may be special cases were 
the prefactors a or b are suppressed by the specific choice 
of the unit cell. The triangulated network is one example 
where a = and the bending mode is suppressed. Below 
we will encounter another example when studying the 
square lattice. For these systems the special geometry 
of the unit cell, or more generally, the local architecture 
has to be taken into account. This is indeed the main 
focus of this article. On the other hand, by assuming 
affine displacements no cooperativity between the elastic 
responses of neighboring cells is possible. The macro- 
scopic modulus G directly reflects the elastic properties 
of the single cell. The local geometry is being hidden 
in the prefactors a and b, while the effect of the assem- 
bled structure may simply be predicted by counting the 
numbers of cells. 



FIG. 4: Shear modulus GPj k as a function of l p /l s for a 2d 
honeycomb foam structure with varying degree of randomness 
A. The blue curve ("Random") corresponds to a "maximally" 
random foam generated from a Poisson point process. Inset: 
Distribution P of segment lengths for the same systems. At 
low levels of randomness (A = 0.3) it can be approximated 
by a Gaussian probability distribution (dashed line), while it 
shows significant broadening upon increasing the randomness 
to A = 0.6, 0.9. The peak disappears completely in the case 
of the maximally random Poisson foam. 



can be inferred from Eq.([9]). By defining 

h * Kip/% , (10) 

and substituting kj 1 — > k^ 1 + on e finds for the mod- 
ulus 



a + b[ — s - 



h 

T 



(11) 



where we have inserted Eqs.([8|) and (|10|) and used the 
relation R — l p /r. This has to be compared with the 
actual results of our numerical analysis in FigQJ The 
normalized shear modulus Gl^/n is shown as a function 
of persistence length l p /l s expressed in units of the av- 
erage segment length. The curves correspond to varying 
degrees of randomness A. 

We find that regular networks (black curve, circles) 
characterized by a single mesh-size / s indeed display the 
functional form expressed through Eq. ljlip . For mesh- 
sizes much larger than the persistence length l s l p 
the network deforms by pulling out thermal undulations 
and G oc k± oc lj A (left part of Fig0|). Decreasing the 
mesh-size beyond l s ~ r stretching of the polymer back- 
bone dominates the modulus G oc k s oc I J 1 (right part 
of Figf?]) . The physically relevant situation for studying 
stiff polymers, however, corresponds to the intermedi- 
ate regime, where the persistence length is much larger 
than the mesh-size, which is still much larger than the 
polymer radius l p 3> l s S> r. Typical actin networks 
modes with stiffnesses k s ,kn and respectively. The with l p = 17/im and r = 5nm may have mesh-sizes in 



C. Cell Polydispersity 

We have tested the validity of the affine model in a 
simple two-dimensional cellular structure with varying 
degree of randomness. We have taken the seeds for 
a Voronoi construction of a regular, honeycomb lattice 
structure and randomly displaced them with a uniform 
probability distribution of width A • l s . The influence of 
randomness on the elastic properties of mechanical (non- 
fluctuating) foams has been studied extensively by var- 
ious authors [H, [H, [HI]. Here, we also include effects 
from thermal fluctuations such that the response of a 
polymer segment is characterized by three deformation 
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affine prediction for the modulus of this system (d = 2) the sub-micron range l s « lOOnm. In this regime, most 
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of the energy is stored in the bending modes leading to 
G oc k± oc l~ 3 corresponding to the plateau region visible 
in FigHl 

Using the values a = 0.2, b = 0.35 and c = 0.14 we 
managed to fit the scaling function of Eq. (fTTj) to the nu- 
merical data (in fact, this is the dashed line in FigJ3]). 
Increasing the level of randomness the presence of the 
additional variable A spoils the scaling property and a 
fit is no longer possible. The power law regimes gradually 
shrink and the cross-over regions increase in size. While 
the mechanical stretching regime is hardly affected by the 
randomness at all, this effect is most pronounced in the 
cross-over from the bending to the thermal stretching 
dominated regime. The physically most relevant inter- 
mediate plateau regime disappears completely and shows 
strong amplitude modulations. 

We have also generated foams by Voronoi tessellation 
of a fully random distribution of points, corresponding to 
a Poisson process (blue curve, left triangles). For these 
"maximally random foams" one could rather use an ex- 

- — 7/2 

pression G oc l s to characterize the modulus at these 
intermediate parameter values. At this point this is only 
an empirical observation. Later, in the context of the fi- 
brous architecture, we will see how this exponent can be 
derived from a scaling argument that properly takes into 
account the randomness in the system. 

One may infer from the inset of Figf?] that deviations 
from the scaling form presented in Eq. pip are indeed in- 
timately connected to a broadening of the segment length 
distribution P(l s ). In the regular structure the distri- 
bution can very well be described by a Gaussian cen- 
tered around the average mesh-size l s (dashed line in 
FigH]) . Random foams, on the contrary, display signifi- 
cantly broader distributions and even have non-negligible 
weight on very small segments. 

We will see below that the different effect of random- 
ness in the thermal and the mechanical stretching regimes 
can be traced back to the unusually strong length de- 
pendence of the entropic stretching stiffness fen oc Zj 4 as 
compared to k s oc Ij 1 . We will find that this leads to the 
breakdown of the affine model whenever there is a suffi- 
ciently broad distribution of segment lengths. Thermal 
networks are thus inherently more sensitive to elements 
of randomness than purely mechanical systems. 

It is instructive to consider yet another lattice struc- 
ture as a basis for our foam model (see Fig.[TJ) . By placing 
the Voronoi points on a slightly randomized square lattice 
one can generate a foam with a bimodal segment length 
distribution having a second peak at some small length 
li (see inset Fig||)]). To understand this, one has to re- 
alize that a generic foam structure generated by Voronoi 
tessellation has only three-fold connected vertices, while 
they are four-fold connected in the square network. A 
small amount of randomness therefore induces a bifurca- 
tion of a four-fold vertex into a short segment with three- 
fold connected vertices at its ends (see FigJ5J). Unlike 
the honeycomb foam, the resulting structure is elastically 
anisotropic and has 3 distinct moduli 36] . In addition to 
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FIG. 5: (a) Pure shear deformation of the square lattice and 
(b) simple shear. Illustration of the bifurcation leading from 
a four-fold connected vertex to a three-fold connected one by 
introducing small amounts of randomness Al a . 
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FIG. 6: The two different shear moduli for the slightly ran- 
domized square lattice (A = 0.3) as shown in Fig[TJ Also 
shown is the modulus of the highly random Poisson foam. 
Inset: Distribution of segment lengths for the same network. 
The dashed line is a fit to a sum of two Gaussians centered 
around h/l s = 0.092 and h/h = 1-35. 



the bulk modulus there are two shear moduli correspond- 
ing to simple and pure shear deformations. These two 
modes are schematized in FigJSJ while the correspond- 
ing moduli (together with the isotropic shear modulus 
of the Poisson foam) are shown in Figl6] Pure shear 
leads to deformations along the main axis of symmetry 
of the unit squares and thus to stretching of the elements. 
The bending regime is therefore strongly suppressed. On 
the other hand, simple shear deforms the squares along 
their diagonals and thus favors the bending mode. Only 
when the stretching energy stored in the small segments 
wh — fcii (h)S a g (h) 2 oc lpH\ equals the bending energy 
in the average segment w± = kj_(l s )S a ff(l s ) 2 oc Ij 1 does 
the system cross-over to a stretching dominated network. 
Noting (from the inset of Fig[6]) that h rj Z s /10 we 
find that this happens when l p k, 10~ 2 Z s in accord with 
FiglU It is interesting to see that the network looses its 
anisotropy at the two points l s = L and l s — r, where 
the modulus takes the same value as that of the Poisson 
foam. This follows from the fact that the stiffness of the 
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average polymer segments is isotropic at these parame- 
ter values and either ku w k± or k s w k±. Comparing 
absolute values we find that the shear modulus in the 
thermal regime, strongly influenced by the presence of 
the small segments, can vary orders of magnitude while 
the mechanical stretching regime is hardly affected at all. 

To conclude this section we emphasize once again that 
polydispersity in the segment lengths can have strong ef- 
fects on the macroscopic elasticity of a cellular polymer 
network. It can lead to modifications of the scaling prop- 
erties, as we have found in the most random foams, as 
well as to quantitative changes of the modulus by several 
orders of magnitude as in the anisotropic square struc- 
ture. As a consequence, experiments which are limited to 
restricted parameter windows would most likely measure 
effective exponents that lie in between the extremal val- 
ues given by pure stretching and bending. One, therefore, 
has to be cautious interpreting experimental data within 
the context of the foam-model, without the knowledge of 
the polydispersity of the structure. 

IV. FIBROUS ARCHITECTURE 

Looking at pictures of cross-linked actin networks re- 
constituted in vitro [l2l [l5| one might wonder whether a 
description in terms of a cellular architecture is actually 
relevant for these systems at all. Besides having a strong 
polydispersity in cell sizes, real polymer networks seem 
to have a hierarchical architecture that allows for smaller 
cells to be generated within larger cells within even larger 
cells. On the contrary, foams only have one of these hi- 
erarchies (see FigJTJ) . What is more, cellular structures 
do not account for the effects of the polymer length If, 
which constitutes an additional mesoscopic scale in the 
problem. 

In the following we want to quantify the effects of 
the polydispersity in connection with the length scale 
If by studying the elastic properties of a generic two- 
dimensional fibrous structure which is defined as follows. 
N anisotropic elastic elements, geometrically represented 
by straight lines of length I f , are placed on a plane of area 
A = L 2 such that both position and orientation of the 
elements are uniformly random distributed. This ran- 
domness entails a distribution of angles 9 e [0, 7r] between 
two intersecting filaments 

P(0) = ^, (12) 

which has a maximum for filaments at right angles. At 
any intersection a permanent cross-link with zero exten- 
sibility is generated. This constrains the relative trans- 
lational motion of the two filaments. For the rotational 
degree of freedom one may introduce an energy contribu- 
tion W IQ t = m{4> — 0o ) 2 for the change of relative cross- 
link angles <p from their initial values (f>o- We restrict 
ourselves to the study of the two limiting cases, where 
the potential is either soft (m — » 0) and therefore allows 



for free relative rotations of the filaments (free hinges), 
or infinitely stiff (m — > oo) and inhibits any change of 
the angles at the cross-links (fixed angles). 

The remaining elastic building blocks of the network, 
the polymer segments, span the distance between two 
neighboring cross-links on the same polymer. Their 
length can be shown to follow an exponential distribu- 
tion 

P(l s ) = l^e- 1 */ 1 * . (13) 

The mean value l s is given in terms of the density p — 
Nl/A as 

T.=*/2p, (14) 

which is a realization of Eq.©. On average there are, 
thus, x = If /l s ~ Ifp segments per polymer. The sim- 
plicity of this network, which has only one structural pa- 
rameter p, makes it an ideal candidate to obtain physical 
insight into the relation between architecture and elastic 
properties of the constituents. This model has frequently 
been used to study the elastic and brittle properties of 
athermal paper sheets [25|, [H, H^, |4(j. In the context 
of biological networks of stiff polymers it has been in- 
troduced in [24[ and recently studied by various authors 
0,113, In ah this work, however, the elastic proper- 
ties of the polymers are modeled by the classical theory 
of Euler-Bernoulli beams. Here, we concentrate on the 
effects of thermal fluctuations, a brief account of which 
we have published recently [231 ]. 



A. Simulation Results 

In Figs.[JJand[5]the results of our simulations are shown 
for fibrous networks with a varying number x of cross- 
links per polymer. The axes are the same as in previous 
plots. The normalized shear modulus Gl^ /k is shown as a 
function of persistence length y = l p /l s expressed in units 
of the average segment length. Short fibers with few 
cross-links, corresponding to low densities are depicted 
in Fig. [7J long fibers or high densities in Fig. [5) In both 
figures we find a regime at large values of the persistence 
length lp/l s (right part of the plot) where the dimension- 
less shear modulus decreases as G oc l~ 2 oc r~ 2 . This cor- 
responds to a purely mechanical stretching regime where 
G (xk s consistent with the mean-field picture of Eq. (TTT|) 

0, S3, Hi m . 

Our main interest, however, lies in the regime of 
lp/L < 10 3 , where the persistence length is small enough 
for thermal fluctuations to become relevant. In this 
regime one may safely neglect the mechanical stretching 
stiffness and set k s — > oo. Then, dimensional analysis for 
the shear modulus requires 

G=frg(x,y), (15) 
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FIG. 7: Scaling function g as a function of y — l p /l s for var- 
ious x — If /l s < 30 (from networks with freely hinged cross- 
links); For rather small values of x = 8, 9 the curves resemble 
the results from the cellular networks. At intermediate values 
10 -3 < y < 10 2 the modulus shows strong modulation and 
develops a dip with increasing x. 
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FIG. 8: Scaling function g as a function of y — l p /l s for 
various x = lf/l s > 60; In contrast to FigQno dependence 
on x is observed any more and new scaling regimes emerge. 
The two branches in the cross-link dominated regime (y <C 1) 
correspond to freely hinged (CLf rco ) or fixed (CLfi xe d) cross- 
link angles, respectively. For i/> 1 one encounters a universal 
thermal regime (U) independent of the cross-link properties 
as well as a mechanical regime (M). 



where we have introduced the scaling variables 

x = lf/l s ~lfp, y = lp/l s ~l p p. (16) 

Comparing to Eq.Q there is an additional argument in 
the scaling function g, the polymer length x = lf/l s - 
This purely geometrical variable counts the number of 
cross-links (or cquivalcntly: segments) per filament. The 
second argument may be written in the alternative form 
y ~ k\\/k±. It characterizes the relative stiffness of 
stretching and bending mode for a typical segment. 
From Fig.[7]one infers that for low densities g = yf(x), 



implying for the modulus G = kuf(plf). This linear de- 
pendence on the "pre-averaged" stretching compliance k\\ 
hints at an entropic stretching dominated regime similar 
to that found in the cellular structures discussed above. 
This regime has been suggested in [2l|, where a scaling 
argument is developed relying on the affine assumption 
borrowed from the mechanical stretching regime. Our 
analysis shows that the domain of validity of this linear 
regime is extremely narrow and confined to short fila- 
ments x < 20 and persistence lengths y <C 1. What 
is more, due to the non-trivial density-dependence ex- 
pressed through the function f(x), the modulus does not 
even display a power-law behavior in the density. In- 
stead, we find that the modulus shows complex depen- 
dence on its variables and develops a dip in the inter- 
mediate parameter region where 10 -3 < y < 10 2 . This 
is also the relevant parameter range for networks of F- 
actin, where the ratio of persistence length to mesh-size 
lp/l s w 10 - 100. 

For medium and high densities Fig.[8]shows non-trivial 
scaling regimes where the scaling function g becomes in- 
dependent of x and therefore of the filament length If. 
This highly non-trivial observation has important impli- 
cations and allows the system to exhibit power law be- 
havior g oc y z . We find non-trivial fractional exponents 
z = 0.46(0.07) and z = 0.9 for small and large values 
of y, respectively. In the figure one can distinguish four 
branches that belong to different realizations of the net- 
work. While branch M (mechanical regime, G oc k s ) has 
been discussed already, the remaining three are obtained 
by setting k s — > oo. The two branches found at small val- 
ues y C 1 relate to networks where the cross-link angles 
are either free to vary (CLf rco , z = 0.46) or are perfectly 
fixed to their initial values (CLfi xc d, z = 0.07), respec- 
tively. We term this regime "cross-link dominated" since 
tuning the cross-link properties may have strong effects 
on the elastic modulus by driving the system from one 
branch towards the other. Both branches merge at y « 1 
where we enter a universal regime (branch U, z = 0.9) 
which is completely independent of the elasticity of the 
cross-links and which therefore is termed "filament dom- 
inated" . 

In all cases, the modulus can be written as a general- 
ized geometric average 

G oc k x - z k\ , (17) 

which has to be contrasted with Eq. (fTTj) . where bending 
and stretching modes are assumed to superimpose lin- 
early (see Table U for a direct comparison of the various 
regimes). There, the system is described either by z — 0, 
if bending dominates, or by z = 1 if stretching is the main 
deformation mode. Values different from the two limit- 
ing cases z — 0, 1 cannot be described by the mean-field 
approach, hence the assumption of affine deformations 
applied on the level of the polymer segments (or the cell 
size) necessarily has to fail. This will become especially 
clear in the following section, where we review the appli- 
cation of affine theories to fibrous architectures. We will 
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TABLE I: Compilation of the different elastic regimes of the 
fibrous network. The modulus is given by G ~ kj^ z kf\ with 
the appropriate values for the exponent z. For comparison 
also the predictions from the theoretical analysis (see below) 
as well as the exponents for the foam structure are given. 
The latter only for fixed cross- link angles (CLfl xc( j), which 
is necessary to make the structure stable. The mechanical 
regime M corresponds to the exponent 2 = 1, however with 
fcii substituted by k s . 



illustrate its failure and highlight the physical principles 
involved. To go beyond we will introduce a model that 
accounts for the spatial distribution of cross-links along 
the backbone of a typical polymer filament instead of just 
considering a single typical polymer segment. This new 
approach will allow us to understand all the features of 
the macroscopic elasticity encountered in Fig[8l 



B. Affine Models in Fibrous Architectures 

In some of the earlier approaches to describe the elastic 
moduli of stiff polymer networks the assumption of affine 
deformations has been applied on the level of the average 
segment which can be characterized by "pre-averaged" 
response coefficients (k(l s )) — > k = k(l s ) introduced in 
Eqs.© and fTU|). The characteristic fibrous structure of 
stiff polymer networks is not accounted for and effectively 
substituted by a highly regular cellular structure. The 
modulus in the thermal regime is then obtained simply 
by replacing in Eq. §§§ the mechanical stretching response 
k s with its thermal counterpart fcii . Several variants of 
this model have been considered in the literature [13, 
HH, |4l| that only differ in the specific (ad hoc) choice of 
the prefactors a, b. The stretching dominated model (28| 
(setting a = in Eq.®) with a modulus depending on 
density as 



G|| ~ p 



(2+d)/(d-l) 



(18) 



and its extensions to nonlinear elasticity [42| , have widely 
been used to fit experimental data for the plateau modu- 
lus in cross- linked F-actin networks [11, 0, 5H • Despite 
this apparent success, it is not clear a priori why in the 
parameter regime of interest the mesh-work should de- 
form by the stretching of bonds when actually bending is 
by far the softer mode — lp/T s 3> 1). In general, 

such a regime can only occur if the specific architecture 
suppresses the soft bending modes as in the triangulated 
structure with its highly coordinated vertices. A second 
approach seems to repair this deficiency by setting in 
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FIG. 9: Illustration of the effects of non-zero residual forces. 
The relaxation of the small segment h from its overly stressed 
state goes to the cost of additional deformations in its neigh- 
bors. Depending on the value of the persistence length the en- 
ergy will mainly be stored in a) the stretching mode (l p ^> l s ) 
or b) the bending mode (l p <C L). 



Eq.© b — 0. The modulus in this theory 



(i+d)/(d-i) 



(19) 



only differs by a factor of p 1 /^- 1 ) from the stretching 
dominated modulus of Eq. (|18p . However, neither the- 
ory provides justification for neglecting the effects of the 
polydispersity in the fibrous system. In fact, if one ex- 
tends the approach to include the distribution of segment 
lengths, such theories necessarily have to fail, as we will 
explain in the next section. 



C. Effects of the Segment Length Distribution 

To understand the origin of this failure consider an 
affine deformation field S^g oc jl s being applied to a ran- 
dom network of stiff polymers with a distribution P(l s ) 
of segment lengths l s . The axial forces f\\ generated by 
such a deformation field can simply be calculated by mul- 
tiplying the deformation with the stretching stiffness of 
the segment 



f\\ = k \\ 5 > 



aff 



Kl pl /ll 



(20) 



Note that in contrast to the purely mechanical situation, 
where the axial force f s = fc s (W ~ nj/r 2 is indepen- 
dent of length, /|| strongly increases with shortening the 
segment length. This implies that, in general, two neigh- 
boring segments on the same filament produce a net force 
Sf at their common node that has to be taken up by the 
crossing filament. There, it leads to additional defor- 
mations that eventually destroy the affine order. This 
mechanism is illustrated in FigJ5] where the relaxation of 
the small segment l\ leads to bending of its neighbor on 
the crossing filament (b). Also the segment I2 on the 
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same filament is affected by the relaxation and experi- 
ences an additional stretching contribution (a). Whether 
the available energy is stored in the stretching or the 
bending mode depends crucially on the value of the per- 
sistence length, as is indicated in the figure. 

One may calculate the probability distribution Q for 
residual forces by summing over all segment lengths that 
are consistent with a given force 5f, 

Q(5f) = (S(\f ]l (l 1 )-f ll (l 2 )\-5f)) . (21) 

The averaging procedure defined by the angular brackets 

(A)= J dhj dl 2 P{h,l 2 )A(h,l 2 ), (22) 

involves the two-point probability P(h,h) of finding 
neighboring segments with lengths l\ and respectively. 
In the special case of the random network considered 
here, there are no correlations between neighboring seg- 
ment lengths such that the distribution factorizes. The 
formula can thus be evaluated by substituting /|| oc Z~ 3 
taken from Eo. (f2"0")) . This inverse relationship between 
forces and segment lengths translates the weight of the 
probability distribution Pq = P(l s — > 0) ^ at small 
segment lengths into polynomial (fat) tails of the corre- 
sponding distribution of residual forces 

Q(Sf -► oo) ex 5f- 4/3 Po , (23) 

which has a diverging mean value. The exponent can 
readily be derived from evaluating the integral measure 
df\\ cx l~ i dl oc /,| dl. As a consequence there are always 
residual forces high enough to cause additional deforma- 
tion of the crossing filament. Hence we conclude that an 
affine deformation field is unstable and that the system 
can easily lower its energy by redistributing the stresses 
to relieve shorter segments and remove the tails of the 
residual force distribution Q(Sf). 

Even though we have evaluated Eq. (f2"T]) for the spe- 
cial case of an exponential segment length distribution 
Eg . (flU]) , it is important to note, that the observed sensi- 
tivity is not a special feature of the fibrous architecture, 
but applies to any polymer network with a broad dis- 
tribution of segment lengths independent of the dimen- 
sionality of the network. Due to the strong length depen- 
dence of A; 1 1 (l s ) the thermal response is highly sensitive to 
even small polydispersity as we have already seen in the 
random cellular network of Sect lHI Cl On the contrary, 
these effects are completely absent in purely mechanical 
models and also in models of flexible polymers, where the 
distribution Q(Sf) degenerates into a delta-function peak 
at the value Sf — 0, and explains the robustness of these 
regimes to randomness. 

If we want to include the effects of randomness into a 
microscopic theory we cannot naively apply the conven- 
tional picture of affine deformations on the scale of the 
single segment. This can safely be done only in highly or- 
dered structures like regular cellular materials. Instead, 



we have to adopt a description of the deformations (at 
least) on the larger scale of the complete polymer. In the 
following we therefore consider a typical polymer filament, 
which is composed of a sequence of segments drawn from 
the distribution P(l s )- To describe the elastic properties 
correctly, we will also have to consider the connections 
of the polymer to the surrounding network matrix, in 
addition to the elastic properties of the segments them- 
selves. We may now employ this picture to explain the 
intricate scaling properties of the polymer network in all 
the parameter regimes displayed in FigJS] 



D. Cross-link dominated Regime 

1. Freely hinged cross-links 

We start with the description of the system in the pa- 
rameter region y <C 1 (l p <C l s ), where the properties of 
the cross-links strongly influence the system's response. 
The idea is to impose a virtual affine deformation on ev- 
ery segment and calculate, as a perturbative correction, 
the contribution to the elastic energy resulting from the 
relaxation out of this reference state. This procedure 
will lead to good predictions only when the corrections 
are small and the affine deformations are only weakly 
perturbed. As we will see below, this is the case in the 
parameter region j « 1. However, it will also become 
clear, that a small perturbation for the deformations is 
sufficient to generate completely different scaling prop- 
erties for the macroscopic modulus. For the moment we 
restrict our attention to free relative cross-link rotations 
(branch CLfi xe d), since then the affine reference state is 
particularly simple and contains stretching contributions 
only. 

As explained above any deviation from the affine ref- 
erence state, induced by relaxation of non-zero residual 
forces, will lead to additional deformations in the cross- 
ing filaments. Since it is more likely that two filaments 
cross each other at an angle close to 90°, the induced non- 
affine deformations will mainly be oriented transverse to 
the contour of the crossing filament and are therefore of 
bending character. The value of the exponent z = 0.46 
supports this assumption and indicates that bending and 
stretching deformations in this regime contribute equally 
to the elastic energy even though the bending mode is 
very stiff (k\\/k± ~ y <C 1). Therefore any relaxation 
of residual stretching forces, will be punished by high 
amounts of bending energy (see Fig.[9j3). Only the small- 
est segments on the polymer, corresponding to the out- 
ermost tails of the residual force distribution, will have 
sufficient energy to perturb the deformation field and re- 
lax to a state of lower strain. 

In the following, we will assume that segments up to 
a critical length l c - to be determined self-consistently 
- fully relax from their affine reference state to give all 
their energy to the neighboring segment on the crossing 
filament. The total energy of the polymer can then be 
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calculated from segments with l s > l c only. There are 
two contributions. First, a stretching energy 



w s {l s ) 



4 



(24) 



from the imposed affine strain field <5 a g oc 7/3. Second, a 
bending energy that is due to the relaxation of a neigh- 
boring segment on the crossing filament out of its affine 
reference state. This process requires that the segment 
of length l s moves the distance 8 a R = 7/5, which cor- 
responds to its own affine deformation. The resulting 
bending energy 



w b {l s ) 



I 2 

2 L s 

K1 H 



(25) 



now depends on the length l s of the segment under con- 
sideration as well as on the length l s of the neighbor- 
ing (now relaxed) segment. As we have assumed above, 
the second contribution w& only arises if the length l s 
is shorter than the critical length l c . The total defor- 
mation energy along the polymer is then obtained by 
adding both contributions and integrating over all seg- 
ments Z s > l c along the filament as well as averaging over 
neighbors with l a <l c , 



Wpoi ~ (//p)«7 5 



dl s P{l s )[ l f 2 




where the prefactor Ijp just counts the number of seg- 
ments per polymer. For simplicity, we have not consid- 
ered any dependence of the deformations on the orienta- 
tion relative to the macroscopic strain field. In essence, 
this would only introduce some additional numerical pref- 
actors that are irrelevant for the scaling picture developed 
here. The integrations are reparametrized by introducing 
the non-dimensional variable A = pl s such that we arrive 
at the expression for the average polymer energy 



Wpoi ~ K'y 2 lfp 2 (pl p /X c + A c ) . 



(27) 



where numerical constants have been dropped and A c :— 
pl c C 1 in the parameter range of interest. Minimizing 
with respect to A c determines a new non-affinity length 



pom = ynm. 



(28) 



that sets the maximal scale up to which the destruction 
of affine deformations lead to a lowering of the elastic 
energy. Inserting this length into Eq. (|2~T|) and multiplying 
by the number-density of filaments p/lf one arrives at 
an expression for the modulus G ~ W™™ ■ p/lf J 2 — 

Kp 7 / 2 ^ 2 . Rewriting the result as 



G 



k 1 fcii 



oc 



r-7/2 



(29) 



we immediately see that our theory reproduces the em- 
pirical result of Eq. (|17p with an exponent z = 1/2, which 
compares well with the measured value of z = 0.46. 
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FIG. 10: Fraction of energy stored in the various segment 
lengths; the curves correspond to different persistence lengths 
at a density of lf/l s = 80, equivalent to l s /l ~ 2 ■ 10 -2 . 



The non-trivial behavior of G observed in Fig. [5] can 
thus be explained by a non-affinity length scale l™ m ~ 

■\Jl s lp below which the affinity of the deformation field 
breaks down. Recapitulating the results from the cellular 

networks in Fig[4j we observe that the same intermediate 

— 7/2 

scaling behavior of G oc l s is found in both architec- 
tures. We have thus established the microscopic origin 
of the scaling law. It derives from a continuous unload- 
ing of smaller segments driven by an interplay between 
segment length distribution and elastic properties of the 
single polymer. This mechanism is illustrated in Fig. ITTJ1 
where a histogram for the fraction of energy stored in 
segments of various lengths is shown. For very small 
persistence length, a significant fraction of the energy is 
stored in the shortest segments. Affine deformations can 
be seen as a good approximation. Increasing the persis- 
tence length, the short segments one after the other loose 
their energies in favor of additional excitations in longer 
segments. This is fully consistent with the assumption of 
a growing non-affinity scale /™ m below which no energy 
is stored. 

It is important to realize that our derivation of the 
exponent does not make use of the precise form of the 
segment length distribution P{l s )- In fact, there is no 
need to perform the integrations explicitly and only the 
limiting behavior of P(l s — > 0) enters. Thus, the conclu- 
sions are valid for a general class of functions that may 
even be slowly vanishing at zero segment length. 

We have also conducted simulations that assume a 
more general form for the stretching stiffness 



1° 

fc llO) = 6k 73T^ 



(30) 



which reduces to the original definition for a = 1. Since 
the relative stiffness of the deformation modes is now 
kn/k± oc (l p /l s ) a we can think of the phenomenological 
exponent a to tune the anisotropy of the individual seg- 
ment. It allows us to extend our discussion to the broad 



12 



10" f 



^ 10 



a 1 


1 1 L<5 
/O ; 







10" 




FIG. 11: Scaling function g as a function of y = Z p /Z s ; the 
symbols correspond to values of a — 0.5, 1.5. In addition the 
scaling function for a = 1 taken from Fig(8]is shown (dashed 
line). Inset: Exponents z determined from the slopes of the 
branches CLf ree (squares) and U (circles), respectively. The 
solid lines correspond to the curves z — a and z = a/(l + a) 
as derived in the main text. 



class of systems for which fcy is a monomial (with units 
energy per area) involving one additional material length 
l p . Repeating the scaling theory for general values of a 
gives z (a) = a/(X + a) which is verified by the results of 
the simulations presented in Fig llll It provides further 
evidence for the validity of our scaling picture. 



2. Fixed Cross-link Angles 

If we want to apply the same reasoning to the network 
with the fixed cross-link angles, we face the problem that 
even a perfectly affine displacement of all the cross-links 
induces some amount of bending of the segments, in ad- 
dition to the usual contribution from the stretching de- 
formations. While an affine strain 7 would change all 
angles by an amount Aip oc 7, due to the infinite ro- 
tational stiffness in the cross-links this cannot actually 
occur. The segments therefore have to experience an ex- 
tra bending contribution induced by cross-link rotations 
(—A(j>) that restore the angles to their original values. In 
the parameter regime y <C 1, where bending is the stiffer 
mode, we therefore expect strong contributions to the 
energy from the bending mode already in the affine ref- 
erence state. Allowing for the relaxation of the smallest 
segments from their stretched state to even stronger re- 
duce the amount of stretching we might find an exponent 
as low as z — 0.07, signalling nearly exclusive contribu- 
tions from the bending mode, not too surprising. In fact, 
we will argue below that neglecting the stretching en- 
ergies, i.e. assuming an exponent of z — 0, represents a 
reasonable approximation to the elasticity in this regime. 



E. Universal Regime 

By increasing y from its small value we soften the bend- 
ing mode and therefore reduce the influence of the fixed 
cross-link angles on the elastic energy. At the same time 
the non-affinity scale Z™ m oc y 1 ^ 2 increases, indicating 
ever stronger deviations from the affine reference state. 



When, eventually, Z" 



I. (K 



y 



1) the 



affine strain field does not serve as a reference config- 
uration any more, since it is strongly perturbed by a 
majority of segments with l s < Z™ m . At this stage, the 
two branches, present in the cross-link dominated regime, 
merge and one enters a universal filament-dominated 
regime. There, the specific properties of the cross-links 
do not influence the macroscopic elasticity notably. 

While the scaling argument presented for the hinged 
network ceases to be applicable, the remaining residual 
forces Sf continue to lead to a redistribution of stresses 
from shorter to longer seg ments, albeit at higher scales. 
As we have shown in 23[, eventually about 90% of the 



energy is stored in the longest 30% of the segments only. 
The new feature as compared to the regime CLf rcc is that 
unloading of a segment from its stretched configuration 
will also lead to stretching of its immediate neighbor on 
the same filament (see Fig (5^). This way, the available 
energy for bending of the crossing filament, which was 
the primary contribution in Eq. I|26p . is reduced. In the 
limit !/ > 1 we can neglect these contributions and cal- 
culate the energy from the polymers' stretching stiffness 
only. The physical picture is that of a serial connection of 
infinitely many segments along the backbone of a "typi- 
cal" polymer. The stiffness of this polymer, and therefore 
the modulus, may be obtained from the stretching spring 
constants of the individual segments ku(l s ) as 



G- 1 = / tfl,P(I,)fc||(J,)~ 1 ocku 



(31) 



corresponding to the exponent z = 1. For the more gen- 
eral response coefficient of Eq. ([30|) this argument predicts 
z — a, a result which is closely confirmed by the results 
of the simulation as can be seen from the inset of FigfTTI 
Note, that the shear modulus in this asymptotic region 
takes the same form as postulated by the affine theory 
in Eq.lfTTj). However, using Eq. ({3"Tj) one can resolve the 
effects according to segment length to find that the con- 
tribution to the total energy from segments with length 
l s grows as W(l s ) oc l*. This strong increase is in accord 
with the assumption of a large non-affinity scale, below 
which no energy is stored, and in striking contrast to the 
affine theory that would yield W a g(l s ) oc fe|| (l s )5^ s oc l~ 2 . 



V. CONCLUSION 

We have studied the macroscopic elastic properties 
of networks of semi-flexible polymers. We provide ex- 
haustive numerical studies supplemented by scaling ar- 
guments that elucidate the subtle interplay between the 
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architecture of the network and the elastic properties of 
its building blocks. 

The main conclusion to be drawn is that, irrespec- 
tive of the specific architecture chosen, thermally fluc- 
tuating stiff polymer networks are inherently more sensi- 
tive to polydispersity and randomness than their purely 
mechanical counterparts. This is due to their strongly 
length-dependent entropic stretching response ku (I) oc 
k 2 // 4 which has to be contrasted with the mechanical 
stretching stiffness k s (l) oc k/1. 

Although simulations have only been conducted in two- 
dimensional networks, the identified mechanism by which 
the structural randomness influences the elastic proper- 
ties is expected to be of universal character and hold 
independent of dimensionality. As we have shown, the 
actual consequences of this susceptibility (e.g. scaling 
behaviour of elastic moduli) may vary from system to 
system and certainly also with the dimension. A precise 
knowledge of the network architecture is therefore indis- 
pensable for the interpretation of experimental data. For 
this it will be most important to develop new techniques 
that allow the characterization of the microstructure and 
monitor its changes upon deformation. As exemplified by 
the discussion in the universal regime Sec. IfVEl where 
the (non-affine) elastic modulus turns out to be similar 
to that in an affine theory, we have shown that macro- 
scopic measurements alone do not suffice to extract the 
network mechanics also on the microscopic scale. 

We have described how the polymer length If can be 
used to drive the system from a simple cellular struc- 
ture with filaments as short as the mesh-size If « l s , to 
a fully scale-invariant fibrous structure characterized by 
infinitely long filaments If — > oo. Especially the latter 
limit allows for intricate scaling behavior that impres- 
sively demonstrates the qualitative difference between 
thermally fluctuating and purely mechanical elastic net- 
works. 

The elasticity of a simple cellular structure may be de- 
scribed by a serial connection of their elementary defor- 
mation modes bending and stretching, respectively. This 
leads to the modulus of Eq. © 

G- 1 = akl 1 + bk^ 1 . (32) 

In this picture, deformations can be drawn from either 
mode and it will be the softer one that dominates the 
modulus. In fibrous networks with fixed cross-link angles 
we have shown that the modes rather act as if they were 
springs connected in parallel. The modulus can then be 
approximated by 

G = ak ± + bk\\ , (33) 

where the prefactors a, b depend weakly on the scaling 
variable y ~ ku/k± ~ l p /l s - The network elasticity is 
therefore always dominated by the stiffer mode, qualita- 
tively similar to a triangulated network, where the spe- 
cific geometry of the unit cell always imposes stretch- 
ing deformations on the system, no matter how soft the 



bending mode actually is. The fibrous architecture ap- 
parently also suppresses the transition into regimes where 
the softer mode is dominant. This conclusion is consis- 
tent with recent simulations on the purely mechanical 
fiber model [l9l [20| . where a transition into a regime 
dominated by soft bending modes (y 3> 1) could only be 
observed at finite values for the filament length If. In- 
creasing the length to asymptotic values If — > oo, as we 
have done here, such a "bending-soft" regime is strongly 
suppressed and eventually cannot occur any more. In- 
stead, the elasticity is governed by the much stiffer (me- 
chanical) stretching mode. A detailed theoretical expla- 
nation of how this suppression is generated in mechanical 
fiber networks will appear elsewhere [44], however, it is 
clear that the mechanism that leads to bending in cel- 
lular structures cannot work in fibrous networks. The 
fact, that any segment is part of the larger structure of 
the polymer fiber leads to strong geometric correlations 
and imposes very strict conditions on possible segmental 
deformations. 

Allowing the filaments to freely rotate at the cross- 
links, a situation which may be relevant for F-actin net- 
works cross-linked for example with a-actinin, we also 
find an asymptotic scaling regime where stretching and 
bending modes contribute equally to the elastic energy, 
Eq.(nj, 

G oc k x ^ z k\ . (34) 

By quantifying the degree of co-operation between neigh- 
boring elements in the network we were able to identify 
a non-affinity length-scale l c below which the state of 
affine deformations is rendered unstable. A scaling ar- 
gument is supplied that allows the calculation of the ef- 
fective macroscopic exponents starting from this micro- 
scopic picture. 

It seems that the effects described above can only be 
accounted for by going beyond the conventional approach 
that considers typical polymer segments only. Instead, 
we propose to describe the elasticity in terms of a typical 
polymer filament and the spatial distribution of cross- 
links along its backbone. By controlling the architecture 
of the network, the scale of the polymer length If there- 
fore seems to implicitly influence the elastic properties of 
the system even in parameter regions where it does not 
enter the macroscopic elastic moduli explicitly. 

Acknowledgments 

We gratefully acknowledge fruitful discussions with 
Mark Bathe, Oskar Hallatscheck and Klaus Kroy. 

APPENDIX A: STIFFNESS MATRIX 

This appendix derives an expression for the stiff- 
ness matrix of a polymer segment imbedded in a two- 
dimensional network. 
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The differential equation governing the bending of a 
beam of length I is given by kX^ — 0, where the trans- 
verse deflection X is induced by the forces Fq, Fi as well 
as the torques Mo, Mi acting on both ends. The solution 
can then be written as 

X(s) = X + X^s + ^-(M -sF /3), (Al) 
Zk 

while equilibrium conditions require that 

F = -F , Mi = - (Mo -F l). (A2) 

Stretching the beam to the position Z is governed by 
the equation 

Z(s) =Z + s- 4-T , (A3) 



with the condition 



Ti = -T Q , (A4) 



balancing the axial forces T. 

The two variables (X, Z) are the coordinates (in the 
frame of the fiber) of the vector u introduced in the 
main text. The rotation is given by 9 — X' . The four 
Eqs. ((AT]) . (fA"2j) , (|A"3"]) and (|AT| can now be inverted to 
yield the forces in terms of the displacements at the beam 
ends (cross-links) 
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where we have defined A = l 2 A/I = 4(//r) 2 The sec- 
ond equality only holds for circular beam cross-sections, 
where the moment of area / = 7rr 4 /4. The corresponding 
matrix is called the stiffness matrix. 

If, in addition to Eq. (|A3|) , we assume that the stretch- 
ing response is governed by that of a thermally fluctu- 
ating stiff polymer we have to take into account fc|| of 
Eq.(|4]). This is achieved by letting both stretching modes 



act in series and substitute k s 1 — > k s 1 + fc|| 1 . Equiva- 
lently, one can assign an effective polymer radius 

4Z 3 

rial = r 2 + — , (A6) 

which now depends on the segment length I as well as on 
the persistence length l p of the polymer. 
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